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Abstract 

We show by means of the theory of the order parameter phase fluctuations that the temperature 
of the "closing" (or "opening") of the gap (and pseudogap) in the electron spectra of supercon- 
ductors with anisotropic order parameter actually takes place within a finite temperature range. 
Every Fourier-component of the order parameter has its own critical temperature. 
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It is well known that the it is not easy to build a self-consistent theory of the copper oxide 
high-temperature superconductors (HTSCs) due to the necessity to take into account differ- 
ent properties of the cuprates, in particular strong electron correlations, low dimensionality 
of the electronic and magnetic properties, anisotropy of the order parameter, pseudogap 
presence, disorder etc. It is extremely difficult to include all these properties into the theory 
self-consistently. The choice of the properties is usually dictated by the aim of studies. Below 
we make an attempt to show that it is possible to explain such an unusual phenomenon as 
smooth disappearing ( "closing" ) of the pseudogap along the Fermi surface arcs from nodal 
points to M-points ((0, 7r) or (tt,0)) of the Brillouin zone in the momentum space. 

It has been proved that the description of the superconductivity in the 2D metals with 
arbitrary carrier density requires not two, as usual, but three self-consistent equations (see 
review Two of them are well known. The first one is the gap equation, which defines 
the order parameter. It can be written in the following form in the case of a separable 
interaction Q]: 

Y 1 r +oc Ap 

A, = ^(k) X>(q) S j_ M) W mTT[fM (1) 

where 

G q K) = , — p — - (2) 

iuj n - f (k)Ti - A k r 3 

is the Matsubara fermion Green's function in the Nambu representation, and 
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f (k) = — — (cos(k x a) + cos{k y a)) - fi (3) 

is a spectrum of one-particle fermi excitations on the square lattice with the constant a and 
with the bandwidth W. The spectrum energy is measured with respect to the chemical 
potential [i. The order parameter in different pairing channels is A(k) = A/(T)7/(k), where 
I stands for the symmetry of the order parameter. In the s-wave pairing channel we choose 
7s(k) = 1, in the p-wave channel 7p (k) = sin(k x a), and, finally, in the related to the HTSCs 
(i-wave pairing channel 7^(k) = cos(ak x ) — cos(ak y ). u n = (2n + l)nT are the Matsubara 
fermi-frequency, fj are the Pauli matrices, and V is the parameter which characterizes the 
fermion attraction (the separable pairing potential is chosen in the following form: VJ(k, q) = 

7 7 i(k)7l(q)). 
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The second equation can be written as: 

i i r+oa j 

i -m^L ^m^ ImTr[fAie)] - (4) 
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This equation allows to connect the number of the mobile (doped) carriers in the conducting 
band with the dispersion law (J3J) and the chemical potential (usually, holes are the carriers 
in the copper oxide HTSCs). Generally, the value of the chemical potential is not equal to 
the value of the Fermi energy (see Q]). 

After the integration over the energy e and by using (j2J) and (jSJ), equations (JTJ) and (HJ) 
can be transformed to a more familiar form: 

1 V f d2k , r ^ tanh[E(k)/(2T)] . 

1 = V J (2^ 7i(k) 2£(k) ' (5) 



J (2tt) 2 I E(k) V 2T J 



(6) 



where E(k) = A/£ 2 (k) + |A(k)| 2 is the excitation energy of the quasi-particles in the su- 
perconducting state. The function E(k) is equal to zero at so called nodal points (with 
\k x \ = \k y \ = \kp\, where is the Fermi momentum). 

The system © and (0) is the self-consistent system of equation of the BCS theory. Its 
solution in the 2D case describes the temperature dependence of the gap (amplitude of the 
order parameter) and of the chemical potential at given value of rif. It is impossible to 
estimate the correct value of the critical temperature from this system, since the critical 
temperature in the 2D case is equal to zero due to the long-wave fluctuations of the phase 

n 

of the order parameter |3J. However, there is another phase transition in the 2D systems 
- the Berezinskii-Kosterlitz-Thouless (BKT) phase transition with the critical temperature 
T = Tbkt a t which the correlation of the order parameter changes its space dependence 
from the exponential (at T > Tbkt) to the power law (at T < Tbkt)- This transition was 
studied most completely in the case of the spin XF-model with the following Hamiltonian 
jj]: Hxy = -f J2 n n(@n — &h) 2 i where J is the exchange coupling and 8 n , 6^ are the phases of 
the spin vectors on the nearest sites n and n. The critical temperature of the BKT transition 
is defined by the following equation: 

Tbkt = \j (7) 
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The order parameter of the superconducting metal is a complex (two-component) func- 
tion, and usually it can be approximated as following (see, for example 0): $(ri,r 2 ) — 
A(r)e i6, ( R \ where r = ri — r 2 and R = (ri + r 2 ) are the relative coordinate and the center 
of mass coordinate of the pair, correspondingly. In the case of the long-wave approximation 
the kinetic part of the thermodynamic potential has the form of Hxy, but in this case the 
superconducting rigidity plays the role of the exchange parameter. This parameter is a 
function of T c = Tbkt, H and A/(T C ). Equation 



7T 



T c =-J(T c ,fi, Ai(T c )) 



(8) 



makes the system of equations (JHJ), © and (jHJ) closed. This system of equations allows us 
to find the gap, the chemical potential and the critical temperature T c . 

The expression for the function J(T C , fi, A;(T C ) for the case with an anisotropic order 
parameter can be found in the complete analogy with the isotropic s-case Q (see also Q). 
It has the following form: 



w. 



DO 



J(T c ^,A l (T c )) = —T c ^ 

n=—oo 



w 
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; (5 - o). 



16T C> / (2tt) 2 cosh 2 (E(k)/(2T c )). 
This expression together with (jHj) defines the equation for the critical temperature of the 2D 
superconducting metal with arbitrary carrier density: 
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T c = —W 
32 



W 



d 2 k 



(9) 



16T c y (2tt) 2 cosh 2 (£(k)/(2T c )) 

The solution of the system of self-consistent equations (jSJ), ©, and © can be found 
analytically in the s-wave case with quadratic dispersion of the quasi-particles, £(k) ~ k 2 . 
In this case A^(0) ~ y 7 ^/, T c ~ nf and \i is always negative when the fermion density nf —>■ 0. 
In the case case of more general dispersion (see, for example Q ), and, moreover, in the case 
of the anisotropic order parameter, the solution can be found only numerically. This solution 
gives us the dependencies of Tc, \i and A;(T C ) on rif. On the other hand, it is possible to find 
the amplitude A^(T) as a function of T from and The solution of these equations 
at Aj(T) = gives us the carrier density dependence of the critical mean-field temperature 
T C MF in the /-pairing channel. This temperature does not correspond to any observable phase 
transition, since there are no phase transitions in the 2D systems, except 5i^T-transition 



as it was mentioned above This transition is a superconducting phase transition in 
the metal with inter-fermion attraction, despite there is no general spontaneous symmetry 
breaking in all the system (f| (for details, see 

The solution of the system of equations (J5]l. (JBJ), and (jUJ) in the s-wave pairing channel 
is presented in Figs.l and 2. It is important to note that the system of equations can be 
analyzed analytically in this case. It is easy to see the complete symmetry of the solutions 
with respect to the point nf — 1 (the case of half-filling), what was already mentioned in 
0,0]. The two-particle (local) bound states exist at any value of V/W at small values of 
n/ (or 2 — ri/, when we consider the hole pairing). In other words, there is no threshold 
value of the coupling for the bound state formation in the s-wave pairing channel in the 
2D case. In this case the chemical potential is negative at any value of V/W and small 
enough ri/, what indicates the crossover from Bose-Einstein condensation regime to BCS 



superconductivity with carrier density increasing (or coupling decreasing) [1CJ, As it 

follows also from these Figures, the following inequality is always correct T c « T^ IF , and 
the carrier density dependencies are close to those found analytically in the case of the 
quadratic dispersion: T C MF ~ ^/nJ,T c ~ rif in the low carrier density limit. The canonical 
BCS relation 2A S (0)/T C MF = 3.52 holds at any carrier density except very low values of n/ 
(see Fig. 3), but the relation 2A S (0)/T C is increasing with carrier density decreasing, since it 
is easy to see that A s (0)/T c ~ ^Jnj at small nj. 

In the cases of the anisotropic p-wave and <i-wave pairing (Figs. 4-6 and 7-9, correspond- 
ingly) the behavior of the superconducting parameters with doping is qualitatively the same. 
The only difference is that the crossover to superfluidity takes place only above some critical 
values of V . The local pairs with non-zero orbital momentum can be formed only when 



the attraction is strong enough 
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13 |. In the cases of the moderate or low attraction 



the chemical potential of the system is positive at any value of n/, therefore only Cooper 
p- and c?-wave pairs can exist in this case. The chemical potential in this case practically 
coincides with the Fermi energy. As it follows from Figures 5 and 8, the relation T c « T^ F 
holds also in these cases, and the ratios 2A p (0)/T c A/jF and 2Arf(0)/T c MF are even higher in 
comparison with the s-wave case, what shows that the isotropic condensate is more stable 
with respect to the thermal fluctuations. The region between T c and T C MF , where the super- 
conducting fluctuations are incoherent and rapidly decay and the order parameter modulus 
is finite, should be interpreted as the pseudogap region. Above T C MF (which in some papers 
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(for example, in 14, 3|) is called the temperature of the decay of superconducting fluctua- 
tions T sc f) the gap in the superconducting spectrum disappears, while there is no any phase 
transition at the point T C MF or around it. 

There is another principal difference between T c and T C MF due to the anisotropy of the 
electron spectra in the superconducting phase. The temperature T c is critical temperature of 
the phase transition, and therefore it is unique for all the system and for all its excitations. 
This temperature can be measured as the temperature below which the resistivity goes to 
zero, or by the Meissner effect, for example. On the other hand, this statement can not 
be used with respect to T C MF . This temperature is not the critical temperature of any 
phase transition in the 2D system, similarly to the case with the s-wave pairing. The 
important question is whether the critical temperature of the formation of different Fourier- 
components or the order parameter is unique. Even the temperature T c can have such values 
that the following rations are true at different directions of the momentum k: T c < A(k), or 
T c > A(k), including T c 3> A(k). More important question is whether the order parameter 
modulus forms at any value of k. A simple (and, in some sense qualitative) answer to this 
question can be given by using the canonical BCS ratio rewritten in the following form: 

2Az(T = 0) 2A,(T = 0)| 7 *(k)| 

t mf T c MF (k) ' 1 ' 

which shows that every Fourier- component of the anisotropic order parameter has its own 
"closing" critical temperature. It should be emphasized that this temperature does not 
correspond to any phase transition. 

The relation (fTU|) is an estimation which actually defines the temperature T c MF (k) = 
T c MF |7;(k)| (or T sc f(k)). As it follows from this relation, the gaped quasi-particle spectra of 
the 2D superconductor with anisotropic order parameter shows its momentum dependence 
also when the temperature changes. It is important that the gap disappears continuously 
from the nodal point, where 71 (k) = at any T, to the M-points, where the gap and the 
corresponding temperature T C MF are maximal, with temperature increasing. The gapless 
(or, actually "pseudogapless" ) character of the spectra will cover larger parts of the Fermi 
surface with temperature increasing at any values of k, at which T = T / MF (]t). This behavior 

nn 

was observed in the ARPES experiments (see recent review [lfiL Il7| and references therein). 

Thus, it is shown that the self-consistent study of the phase fluctuations of the supercon- 
ducting order parameter allows one to describe qualitatively the anisotropy of the pseudogap 
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and its disappearance within a rather wide range of temperatures. 
Acknowledgments 

One of us (V.M.L.) would like to thank H. Beck, Yu.B. Gaididei and V.P. Gusynin for 
stimulating discussions of the results. The work was partially supported by The Swiss 
National Science Foundation (SCOPES-project 7UKPJ062150.00/1). 



[1] V. Loktev, R. Quick and S. Sharapov, Phys. Rep. 349, 1 (2001). 
[2] J.R. Schrieffer, Theory of Superconductivity, Benjamin, N.Y. (1964). 

[3] A.Z. Patashinskii, V.L. Pokrovsky, Fluctuation Theory of Phase Transitions, Nauka, Moscow 
(1975). 

[4] Yu.A. Izyumov, Yu.N. Skryabin, Statistical Mechanics of Magnetically Ordered Systems, 

Nauka, Moscow (1987). 
[5] V.M. Loktev and V.M. Turkowski, Low Temp. Phys. 30, 179 (2004). 
[6] V.P. Gusynin, V.M. Loktev, and S.G. Sharapov, JETP Lett. 65, 182 (1997). 
[7] It is necessary to mention that usually the temperature T c Mi? is not considered in the case 
of the 2D XY-model. Being rather high (T C MF ~ 4 J for the square lattice) this temperature 
shows when and where (above T^ F ) the spin system can be considered in the paramagnetic 
state, i.e. where the average spin modulus on the lattice is zero. 
[8] E.V. Gorbar, V.P. Gusynin, V.M. Loktev, Fiz. Nizk. Temp. 19, 1171 (1993). 
[9] V.P. Gusynin, V.M. Loktev, and S.G. Sharapov, JETP Lett. 69, 141 (1999). 
[10] M. Randeria, J.-M. Duan and L. Shieh, Phys. Rev. Lett. 62, 981 (1989). 
[11] M. Randeria, in: A. Griffin, D. W. Snoke, and S. Stringari (Eds.), Bose-Einstein Condensation, 

Cambridge U.P, New York, 1995, P.355. 
[12] Yu.B. Gaididei, V.M. Loktev, Ukrainian Journ. of Phys. 19, 253 (1974). 
[13] MYu. Kagan and T.M. Rice, J. Phys. C 6, 3771 (1994). 

[14] J.L. Tallon, in: T. Yamashita and K. Tanabe (Eds.), Advances in Superconductivity, Springer, 
Tokyo, 2000, P. 185. 

[15] S.H. Naqib, J.R. Cooper, J.L. Tallon and C. Panagopoulas, Physica C 387, 365 (2003). 



7 



[16] A. Damascelli, Z. Hussain and Z.X. Shen, Rev. Mod. Phys. 75, 473 (2003). 



[17] M.R. Norman and C. Pepin, preprint cond-mat/0302347 , 



8 




0.5 



-0.5 



b) 


V/W=2.0 - 


1.6 




1.0 




<W ^%; ::? _ 





0.5 



1.5 



n f n f 
FIG. 1: The carrier density dependence of a) A s and b) fi at T = and different values of V in 
the s-wave pairing channel. 
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FIG. 2: The carrier density dependence of a) T^ F and b) T c at different values of V in the s-wave 
pairing channel. 
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FIG. 3: The carrier density dependence of the ratios a) 2A S (T = 0)/T c MF (a)) b) 2A S (T = 0)/T c 

at different values of V in the s-wave pairing channel. 
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FIG. 5: The same as in Fig.2 for the p-wave pairing case. 
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FIG. 6: The same as in Fig.3 for the p-wave pairing case. 
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FIG. 8: The same as in Fig.2 for the ci-wave pairing case. 
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FIG. 9: The carrier density dependence of the ratios a) 4A d (T = 0)/T c MF and b) 4A d (T = 0)/T c 
at different values of V in the d-wave pairing channel (The maximal value of the gap in the ci-wave 
pairing case is 2A^(T = 0), not A^(T = 0) like in the s- and p-wave cases.) 
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